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ABSTRACT 

We present the stellar kinematics of 48 representative elliptical and lenticular galaxies ob- 
tained with our custom-built integral-field spectrograph SAURON operating on the William 
Herschel Telescope. The data were homogeneously processed through a dedicated reduction 
and analysis pipeline. All resulting SAURON datacubes were spatially binned to a constant 
minimum signal-to-noise. We have measured the stellar kinematics with an optimized (penal- 
ized pixel-fitting) routine which fits the spectra in pixel space, via the use of optimal templates, 
and prevents the presence of emission lines to affect the measurements. We have thus gener- 
ated maps of the mean stellar velocity V, the velocity dispersion er, and the Gauss-Hermite 
moments /13 and of the line-of-sight velocity distributions. The maps extend to approxi- 
mately one effective radius. Many objects display kinematic twists, kinematically decoupled 
components, central stellar disks, and other peculiarities, the nature of which will be discussed 
in future papers of this series. 

Key words: galaxies: bulges - galaxies: elliptical and lenticular, cD - galaxies: evolution - 
galaxies: formation - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

We are carrying out a study of the structure of 72 representative 
nearby early-type galaxies and spiral bulges based on measure- 
ments of the two-dimensional (2D) kinematics and line-strengths 
of stars and gas with SAURON, a custom-built panoramic integral- 
field spectrograph for the William Herschel Telescope (WHT), 
La Palma. SAURO N is base d on the TIGER microlens concept 
( Baco n et alJll995t s ee also | Baco n eta! l2000h . and is described 
in detail in Paper I iBacon e^lT]20^^ ). The objectives of th e 
SAURON survey are summarized in Paper II Ide Zeeuw et all2002h . 
which also contains the definition and global properties of the 
sample. Here we present maps of the stellar kinematics for the 
48 elliptical (E) and lenticular (SO) galaxies in the survey. The 



morphology and kinematics of the ionized gas, the line-strength 
distributions, similar measurements for the spiral galaxies in the 
sample, and higher spatial resolution integral-field observations 
obtained with OASIS on the Canada-France-Hawaii Telescope 
(CFHT) will be presented separately. The data and maps presented 
in this paper will be made available via the SAURON WEB page 
http://www.strw.leidenuniv.nl/sauron/ 

In Section [2] we summarize our observational campaign, as 
well as the principles of the data reduction and the measurement 
of the line-of-sight velocity distributions. We present the kinematic 
maps in Section|3| and our conclusions in Sectionfl] 
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Table 1. The SAURON observing runs. 



Run 


Dates 


Clear 


1 


14-20/02/1999 


3.5/7 


2 


08-14/10/1999 


4/8 


3 


27/03-03/04/2000 


3.5/8 


4 


01-04/09/2000 


3/3 


5 


14-25/03/2001 


11.5/12 


6 


14-17/01/2002 


3/4 


7 


10-19/04/2002 


5/10 


8 


24-27/04/2003 


2/4 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observing runs 

We were allocated a total of 56 nights on the WHT, split into 8 runs 
over four years, to observe the complete SAURON representative 
sample of 72 objects. This substantial allocation was divided al- 
most equally between the observing times of the United Kingdom 
and the Netherlands. Table Q summarizes the dates and provides 
weather statistics. A total of 35.5 nights were clear. Run 1 followed 
soon af ter the commissioning period (February 1-5, 1999, see 
IPaper J) . SAURON is a visitor instrument, so each run was preceded 
by one or two days of configuring and calibrating the instrument, 
with help from the Isaac Newton Group staff. The first three runs 
were plagued by poor weather, and technical problems further hin- 
dered progress in Run 4. The installation of ULTRADAS (observing 
control software at the WHT) in the spring of 2001 improved the 
efficiency, which, combined with the exceptional weather statistics 
of Run 5 allowed us to make more rapid progress. A new volume- 
phase holographic grating (VPH) was installed for Run 8, resulting 
in a further increase in the efficiency of the spectrograph to a re- 
markable 20% end-to-end (atmosphere, telescope, detector). 

The SAURON sample of 48 E and SO galaxies is representative 
of nearby bright early-type galaxies (cz ^ 3000 km s ; Mb ^ 
— 18 mag). As described in lPaper III it contains 24 galaxies in each 
of the E and SO subclasses, equally divided between 'cluster' and 
'field' objects (the former defined as belonging to the Virgo cluster, 
the Coma I cloud, and the Leo I group, and the latter being galaxies 
outside these clusters), uniformly covering the plane of elli pticity e 
versus absolute blue magnitude Mb- Tables Al and A2 in lPaper ID 
provide basic information on the objects. 

We used the low resolution mode of SAURON, giving a 
field of view of 33" x 41", fully sampled by 1431 square lenses 
0'.'94 x 0'.'94 in size, each of which produces a spectrum. Another 
146 lenses sample a small region about 1.'9 from the field center 
for simultaneous observation of the sky background. The wave- 
length range 4800-5380 A is covered at 4.2 A spectral resolution 
(FWHM, Uinst = 108 km s" 1 ) with a sampling of 1.1 A per pixel. 
This range includes a number of important stellar absorption lines 
(e.g. H/9, Mg6, Fe) and potential emission lines (e.g. H/3, [Olll], 
[Nl]). 

Each galaxy field was typically exposed for 4 x 1800 s, each 
dithered by a small non-integer number of lenses to avoid system- 
atic errors due to e.g., bad CCD regions. In about one third of the 
cases, we constructed mosaics of two or three pointings to cover 
the galaxy out to about one effective radius R c , or, for the largest 
objects, out to 0.5i? e - The footprints of these pointings are shown 
overlaid on i?-band Digital Sky Survey images in Figure □ For 
each of the 48 galaxies, Table|2|lists the run(s) in which it was ob- 
served, and the exposure times for the individual pointings. Note 



that Run 8 was used to complete the sample of early-type spirals, 
the data of which will be presented in a future paper of this series. 
Arc lamp exposures were taken before and after each 1800s expo- 
sure in order to track flexures and provide a wavelength calibration. 
We also obtained data for a number of velocity, line-strength and 
flux standard stars throughout each night, for calibration purposes. 



2.2 Data reduction 

We reduced the SAURON observations with the dedicated 
XSAURON software developed at CRAL and described in lPaper fl 
The steps include bias and dark subtraction, extraction of the spec- 
tra using a fitted mask model, wavelength calibration, low fre- 
quency flat-fielding, cosmic-ray removal, homogenization of the 
spectral resolution over the field, sky subtraction, and flux calibra- 
tion. The wavelength calibration is accurate to 0.1 A (6 km s _1 ). 



2.3 Merging and binning 

We merged multiple exposures and mosaiced different fields by 
truncating the wavelength domain to a common range, register- 
ing the exposures using reconstructed images, and combining the 
spectra (as well as the noise spectra) with optimal weights and 
(re)normalisation. During this process, we resampled the datacubes 
to a common spatial scale of 0" 8 x 0" 8, reorienting them to have the 
North up and East left , and correcting for the effect of atmospheric 
refraction (as in Emsellem et al. 1996). 

Because a minimum signal-to-noise ratio (S/N) is required 
to reliably measure stellar kinematics, we spatially binned the data 
cubes using an adaptive scheme developed bv lCappellari & Copinl 
(2003). In this approach, the spectra are coadded by starting from 
the highest S/N lenslet, and accreting additional lenslets closest to 
the current bin centroid. A new bin is started each time the target 
S/N is reached. The resulting bin centroids are then used as start- 
ing p oints for a centroidal Voronoi tessellation (see, e.g.. lDu et alJ 
1999), ensuring compact non-overlapping bins and a uniform S/N 
in faint regions. We imposed a minimum S/N of 60 per spec- 
tral resolution element for the kinematical measurements presented 
here. 



2.4 Determination of the point spread function 

We determined the point spread function (PSF) of each merged ex- 
posure by comparing the reconstructed SAURON intensity distribu- 
tion with HST/WFPC2 images. The SAURON PSF is well-described 
by the sum of two circular Gaussians whose parameters are de- 
termined by minimizing the difference between the WPFC2 image 
(convolved by these Gaussians) and the SAURON reconstructed im- 
age fBacon et alj|200l5) . The result is valid in the central region 
where the overlap between individual exposures is maximum. The 
derived seeing values are listed in Table [3] WFPC2/F555W im- 
ages were favoured as the F555W filter corresponds more closely 
to the SAURON spectral domain. WFPC/F814W images were how- 
ever used for about a third of the galaxies presented here, for which 
the F555W filter was not available. HST imaging does not exist for 
two objects in the E/S0 sample, namely NGC 2695 and NGC 4262, 
and for another two, NGC 4387 and NGC 7332, only WFPC1 data 
is available. For these galaxies, we favoured the estimates of the 
PSF width provided by the seeing monitor at WHT. In most cases, 
SAURON undersamples the seeing. 
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Figure 1. i?-band Digital Sky Survey images of all 48 E and SO galaxies in the SAURON representative sample. The size of each image is 4'x4', and 
the orientation is such tha t north is up and east is left. The bar located at the bottom-left corner of each image ind icates the size of 1 R e (from RC3: 
Ide Vaucouleurs et al]|l99ll) for that object where available. NGC 2699, NGC 4477 and NGC 5845 have Re of 9" iBlakeslee etai]|200lh . 26" and 3" 
lAlonso et all2003l) . respectively. Overlaid on each image is the approximate footprint of the SAURON pointings obtained for that object. 
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Table 2. Details on the exposures of the E/SO SAURON representative sample. Characteristics for each target are provided in Paper II <de Zeeuw et all2002h . 
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Notes: (1) NGC number. (2) Run number (see TablelTl. (3) Pointing number. (4) Exposure time, in sec. 



2.5 Stellar kinematics 

Most of the galaxies in this sample are contaminated by sig- 
nificant emission lines of e.g., H/3, the [O III]AA4959,5007 or 
[N I] AA5 198,5200 doublets. This causes problems for techniques 
worki ng in Fourier space, e.g. the Fourier Correlation Quotient 
(FCQ , lBendeJl990l) or the Cross Correlation Fitting (CCF. IStatleil 
1995) techniques, as they do not permit direct masking of specific 
spectral regions. Such methods can obviously make use of continu- 
ous spectral regions not or weakly affected by emission lines. How- 
ever, considering the short spectral domain provided by SAURON, 
and the significant fraction of galaxies found to exhibit the emis- 
sion lines which punctuate this region, these techniques are not 
well suited to these data. For this reason we opted for the use of 
a direct pixel-fitting (PXF) routine to extract the kin ematics, allow- 
ing emission-line regions to be easil y masked (e.g. . Ivan der Marel 
11994 CaDDellari & Emsellem 2004, and references therein). We 
emphasize that all results obtained with this new implementation 
are fully consistent with the early results presented in Paper II. 



2.5.1 Penalized pixel fitting with optimal templates 

We derive the line-of-sight velocity distribution (hereafter LOSVD) 
parametrized by a Ga uss-Hermite function Ivan der Marel & Franxl 
Il993l l Gerhar dl 19931) for each spectrum using the PXF routine: the 
algorithm finds the best fit to a galaxy spectrum (rebinned in In A) 
by convolving a template stellar spectrum with the corresponding 
LOSVD. This provides the mean velocity V and the velocity dis- 
persion a, as well as the higher order Gauss-Hermite moments /13 
and hi, which quantify the asymmetric and symmetric departures 
of the LOSVD from a pure Gaussian (related to the skewness and 
kurtosis respectively). Although this appears as a conceptually sim- 
ple procedure, there are two critical issues (not specific to SAURON) 
to carefully consider to obtain robust and reliable kinematic mea- 
surements. 



Optimal templates: Although robust against the presence of 
emission lines, PXF is sensitive to template mismatch, in contrast 
to CCF and FCQ which minimize the coupling between different 
absorption features by restricting the analysis to the main peak of 
the cross correlation between the template and the galaxy spec- 
tra. The use of PXF thus requires templates which match closely 
the galaxy spectrum under scrutiny. The derivation of optimal stel- 
lar templates must be conducted for every individual spectrum (as 
galaxies may e.g. exhibit significant metallicity or age gradients 
within the SAURON field of view). This is achieved via the use of an 
extensive stellar library spanning a large range of metallicities and 
ages. We primarily included 19 spectr a from the librar y of single- 
metallicity stellar population models o flVazdekisl fl999). and added 
5 spectra from the Jones stellar library (from which Vazdekis' mod- 
els are built; Jonesl l 19971) to provide spectra with strong Mg b in- 
dices. This also requires the inclusion of additive Legendre poly- 
nomials to adapt the continuum shape of the templates. 

Penalized pixel fitting: When measuring the Gauss-Hermite 
moments of the LOSVD up to hi, undersampling becomes im- 
portant when the observed dispersion is less than about 2 pixels 
(equivalent to around 120 km s _1 for SAURON). Below this disper- 
sion a very high S/N is required to accurately measure the Gauss- 
Hermite parameters. At our minimum S/N ~ 60, the data are 
therefore unable to significantly constrain all the (V, a, /13, /14) pa- 
rameters simultaneously in the undersampled regime, and the scat- 
ter in the measurements increases dramatically. In this situation 
one usually wants to reduce the parametric form of the LOSVD 
to a simple Gaussian. We therefore used the penalized pixel-fitting 
(herea fter pPXF) method developed by CaDDellari & E mselleml 
(2004) to perform this biasing of the solution towards a Gaussian 
in a statistically motivated and automatic way, as a function of both 
the S/N of the spectra, and the observed a. 

We implemented both the template optimization and the pe- 
nalization directly into our PXF routine. This works for each indi- 
vidual galaxy spectrum as follows. The direct pixel fitting is a non- 
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Table 3. Characteristics of the E/SO SAURON galaxies. 



NGC Type V syst V bary Source Seeing 9 S 

(I) (2) (3) (4) (5) (6) (7) 



Notes : (1) NGC number. (2) Hubble type (RC3: Ide Vaucouleurs et alj 
Il99ll) (3) Estimate of the heliocentric systemic velocity in kms - 1 (see 
Sect.[3j. (4) Mean barycentric correction applied to the merged datacubes 
(in km s -1 ). (5) Source for the seeing determination (see text) (6) Seeing, 
full width at half maximum in arcsec. (7) Position angle, in degrees, of the 
vertical (upward) axis in the maps shown in Figs. l4al4ll 



linear least-squares problem, for the parameters (V, a, /13, /14) of 
the LOSVD, which is solved using a sequential quadratic program- 
ming method including bounds (NAG routine E04UNF). At each 
iteration we find the best fitting linear combination of the stellar 
templates (with positive weights) convolved with the correspond- 
ing LOSVD. This includes the additive polynomials and provides 
the optimal template for that specific set of kinematic parameters. 
The residual vector r between the observed spectrum and the con- 
volved optimal spectrum is then derived and the penalization term 
is added to obtain a 'perturbed' vector of residuals r', which is fed 
into the nonlinear optimization routine. 

The perturbed vector has the form r' = r + XT)a(r), where 
A is an adjustable parameter of order unity, V 2 = ft| + h\ is the 
integrated squared deviatio n of the LO SVD from its best fitting 
Gaussian Ivan der Marel & Franxlll993l) . and <r(r) is a robust es- 
timate of the standard deviation of the residuals. The qualitative 
interpretation of the perturbation can be illustrated by noting that 
for A = 1, an amplitude of 0.1 for /13 or /14 (10% deviation from 
a Gaussian) requires an improvement of at least 0.5% of the un- 
perturbed <r(r) to be accepted by the optimization routine as a de- 
crease of the perturbed residuals r'. Setting A = provides a non- 
penalized PXF routine. After extensive testing we adopted a value 
of A = 0.7, which makes sure the bias is always small compared 
to the measurem ent errors, when rj > 120 km s - and S/N > 60. 
See lCappellari & Emsellerr] 120041) for details on this pPXF tech- 
nique. 

We estimated errors by a Monte-Carlo method in which the 
kinematic parameters are derived from many realizations of the 
input spectrum obtained by adding Poissonian noise to a model 
galaxy spectrum. Figure [2] summarizes the results from such sim- 
ulations. The significant damping of /13 and /14 towards zero is 
visible when a < 120 km s - , but as noted above, this damp- 
ing effect decreases with increasing S/N and becomes almost 
negligible for S/N > 100, as we have in the center of most 
of our galaxies. For a S/N of 60, the lcr errors on V, a, /13 
and /14 are 4 kms -1 , 5 kms -1 , 0.03 and 0.03 respectively for 
a = 120 kms -1 , and 8 kms -1 , 7 kms -1 , 0.02 and 0.02 re- 
spectively for a = 300 km s . At higher S/N, the errors scale 
roughly inversely proportional to S/N. These estimates are only 
lower limits since they do not account for the effect of template 
and continuum mismatch (Appendix |B3j. Error maps will be made 
available with the public data release. 

2.5.2 The robustness of the kinematic measurements 

In order to assess the robustness of the method described above, we 
conducted an extensive series of tests. We first carefully inspected 
the fits of the observed spectra obtained for a number of galaxies 
and these were always found to be excellent (see Appendix lB2l for 
some illustrations). We checked that the kinematic measurements 
were not affected by the details of the input library. 

Most galaxies required a Legendre polynomial of degree 6 for 
the pixel fitting. For a few galaxies we had to impose a higher de- 
gree (10) to properly fit the spectral shape. We verified that this 
did not artificially bias the kinematic measurements, even the more 
sensitive higher order Gauss-Hermite moments /13 and h^. We also 
found that applying either multiplicative or additive polynomials 
led to perfectly consistent results. 

All tests were performed with two completely independent 
implementations of the pPXF method and both gave indistinguish- 
able results. The stellar kinematics of all 48 galaxies was also ex- 
tracted using FCQ in combination with the gas-cleaning procedure 
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Figure 2. Simulated errors on the measured kinematic parameters using 
Monte Carlo realizations. These plots illustrate the case of an input S/N of 
60, h:j = 0.1, /i4 = —0.1, and for velocity dispersions between 50 and 
360 km s . The solid line in the top panels show the difference between 
the measured (via pPXF) and the input values for the mean velocity V (left) 
and the velocity dispersion cr (right). In the bottom panels, we show the 
recovered /13 (left) and hi (right) values. The dotted lines illustrate the la 
error on these measurements at constant input velocity dispersion. 



described in Paper II. For galaxies not severely affected by gas 
emission the two methods produced fully consistent results (see 
SectionlBTi. 
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Figure 3. Comparison between aperture (i? e /10) measurements of the cen- 
tral dispersion of the SAURON E/S0 galaxies and values taken from the liter- 
ature. Bottom: SA URON centra l disper sions versus valuesfromjDaviey^il 
Il987|) (circles), iFaber et alJ (j.997) (triangles) and Pnmniel & Simien 



1 1996) (squares). Top: difference between the values from SAURON and the 
literature. The errors from both datasets have been summed quadratically. 



2.5.3 Comparison with published kinematics 

We have already shown in Paper II that the SAURON stellar and 
gas kinematics are consistent with published data. To have a more 
global view of the sample, we performed a comparison between 
the central stellar velocity dispersion un from the SAURON data 
and from published values. We made use of three main sources 
which include a significant overlap with the SAURON target list, 
namely the 7 Samurai 1 Davies et alll987l) . the weighted averaged 
values compiled by Faber et al. 1 1997) and the compilation of data 
gathered from papers by Prugniel and Simien and available via the 
Hypercat 1 on-line catalog IPrugniel & SimieiJl 996 ) . Our measure- 
ments were derived by averaging (luminosity-weighted) dispersion 
values within a circular aperture with a radius of 7? e /10, but im- 
posing a minimum aperture of 2.4" x 2.4" (3x3 pixels). A lin- 
ear regression applied to the combined dataset provides a slope of 
0.98±0.03 and an average residual of — 1± 17 km s . The average 
residuals for each individual dataset are all consistent with zero and 
with the scatter expected from the quadratically summed error bars 
(seeFigJ^^OilSkms -1 for the 7 Samurai, -2.1 ± 16 km s" 1 
for Faber et al. (1997) and -5.6 ± 15 km s" 1 for the data of 
Prugniel & Simien ( 1996). There is no measurable trend for sub- 
samples denned either for a < 200 km s _1 or a > 200 km s _1 . 
This comparison demonstrates that our (velocity dispersion) mea- 
surements do not seem to suffer from any detectable systematics. 



1 http://www-obs.univ-lyonl.fr/hypercat/ 



3 OBSERVED KINEMATICS 

Figures I4al41l display maps of the absorption-line kinematics of 
the 48 objects, all plotted on the same angular scale. Not count- 
ing multiple exposures, and taking into account overlap of mo- 
saiced exposures, these maps correspond to a total of about 130000 
galaxy spectra (about 32000 after adaptive binning). The maps are 
displayed according to increasing NGC number. In each case we 
show the total intensity reconstructed from the SAURON spectra, 
the mean stellar velocity V, the velocity dispersion a, as well as 
the Gauss-Hermite moments /13 and /14, as derived from the pPXF 
routine. Mean velocities are with respect to estimated heliocentric 
systemic velocities, the values of which are provided (corrected for 
the barycentric motion) in Table|3| 

The maps display a wealth of structures. Most galaxies in this 
E/S0 sample show a significant amount of rotation. Some galaxies 
exhibit isovelocity contours steeply rising towards the photometric 
major-axis and also tend to show elongated or dumbbell shaped a 
fields, for example NGC 821, NGC 2549, NGC 3377, NGC 3384, 
NGC 4526, NGC 4660, NGC 5308 and NGC 5845. Other objects 
have strongly misaligned photometric and kinematic axes, indica- 
tive of a non-axisymmetric structure (e.g., a bar). This includes 
NGC 474, NGC 2699, NGC 3384, NGC 4262 and NGC 4477. 
Some objects display an increase in a along the photometric major- 
axis associated with a flattening or turnover of the mean veloc- 
ity gradient along the same axis, for example NGC 4550 and 
NGC 4473. Kinematically decoupled components (KDCs), show- 
ing either a twist in the nuclear kinematic axis or a central velocity 
field clearly distinct from the rest of the galaxy, are evident in a 
number of maps, for example NGC 3414, NGC 3608, NGC 4458, 
NGC 5198, NGC 5813, NGC 5831, NGC 5982, and are detected 
at smaller (apparent) scales in NGC 4150, NGC 4382, NGC 4621, 
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Figure 4a. Maps of the stellar kinematics of the 48 E and SO galaxies in the SAURON representati ve sam ple. The SAURON spectra have been spatially binned 
to a minimum S/N of 60 by means of the Voronoi 2D-binning algorithm of Canpellari & Copin 1 2003). All maps are plotted to the same spatial scale. The 
arrow and its associated dash at the top of each column mark the North and East directions, respectively; the corresponding position angle of the vertical 
(upward) axis is provided in Tablel3l From top to bottom: i) reconstructed total intensity, ii) stellar mean velocity V, iii) stellar velocity dispersion <r, iv) and 
v) Gauss-Hermite moments hg and /14. The cuts levels are indicated in a box on the right hand side of each map. 
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Figure 4c. 
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Figure 4f. 
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Figure 4i. 



16 Emsellem et al. 



mC 46£IU-> ]ffGC 4850 ^ NGC 5190 -NGC 5309 




Figure 4j. 



The SA URON project -III 17 




-10 10 i LO i» 



1 8 Emsellem et al. 




Figure 41. 
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NGC 7332 and NGC 7457. Some of these were known previously, 
others are new. Appendix [A] contains a brief description of these 
structures, including references to previous long-slit work on these 
objects. The kinematics presented here will be analysed and dis- 
cussed in detail in subsequent papers of this series. 



4 CONCLUDING REMARKS 

The maps presented in this paper are the result of the first sur- 
vey of the stellar kinematics of nearby early-type galaxies with an 
integral-field spectrograph. This unique dataset demonstrates that 
early-type galaxies display significant and varied structure in their 
kinematic properties, reinforcing the preliminary evidence shown 
in Paper II, and the hints provided by earlier long-slit spectroscopy. 
The present SAURON maps demonstrate that 2D coverage is a pre- 
requisite to properly understand the dynamics of nearby galaxies. 

The maps shown here constitute only part of the information 
contained in the SAURON spectra. The morphology and kinemat- 
ics of the emission-line gas, as well as maps of the various line- 
strengths derived from the same observations will be presented in 
future papers. 
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APPENDIX A: DESCRIPTION FOR INDIVIDUAL 
GALAXIES 

Here we briefly comment on the structures observed in the 
SAURON maps of the E/S0 sample presented in this paper. A 
thorough and quantitative assessment of these structures will be 
carried out in subsequent papers. Independently derived veloc- 
ity and velocity dispersion maps for NGC 3379, NGC 3384, 
N GC 4550 and NGC 7 4 57 w e re presented an d disc ussed 
in lAfa nasiev & Sil'chenko 1 20021). ISil'chenko et alj 120021) and 
ISil'chenko et aljfe003l) . 

NGC 474: This gala xy (Arp 227), fa mous for its shell structures 
(e.g. [Turnbull. Bridge s & CarteJ I 1999), is a clear case of a non- 
axisymmetric object, with a strong misalignment between the kine- 
matic and photometric major-axes <Hau. Balcells & CarteJl996T) . 

NGC 524: This galaxy has surprisingly high velocities (ampli- 
tude larger than 140 kms" 1 ) considering its lo w ellipticity e < 
0.05 jSil'chenkdl2000HSimien & PrugniellE oOO). as well as large 
/13 values up to large radii, anti-correlated with V . 

NGC 821: A close to edge-on galaxy with one of the most promi- 
nent disk-like kinematics in the E/S0 SAURON sample revealed by 
its velocity field (see also |Pinkne vetaljE003l IScorza & Bended 
ll995UBender. Saglia & Gerhard 19941) . 

NGC 1023: This is an example of an SB0 galaxy showing /13 
anti-correlated with V inside ~ 10", and correlated outside 



feoweret all 1200 ll) . There is a strong twist of the zero veloc- 
ity curve nea r the center. Stellar kinemat i cs has been pre viously 
published bv ISimien & Prugniei Jl997bl) . iNeistein et all J 19991) . 
ISil'che nko ( 1999), including multiple long-slit d ata parallel to the 
major-axis iDebattista. Corsini & Aguerril |2002'). The /13 and /14 
maps show significantly negative values on the western side of 
the gal axy which are not observed bv lDebattista. Corsini & Aguerril 
(2002). The cause for this discrepancy is not known, so these fea- 
tures should be confirmed. 

NGC 2549: The central high dispersion region is elongated along 
the photometric minor-axis, due to the contribution of a rapidly ro- 
tating component (see a lso Neist ein et alll999llSimien & Pr ugniel 
Il997bl: ISeifert & Scorzdll996l) . exemplified by the isovelocities 
pinched together along the photometric major-axis, and the corre- 
spondingly large /13 values. 

NGC 268 5: The kinematics of this famous object (the He- 
lix ga laxy; lBurbidge"& Burbidgel ll959t IPeletier & ChristodouloJ 
1993) are strongly perturbed by extinction due to polar dust 
lanes on its North East side. The velocity dispersion in the cen- 
ter displays a double p eak structure 1 Simien &PrugniellH997bt 
lHau. Carter & Balcellsll999l) . 

NGC 2695: This galaxy shows a high amplitude in /13 in the cen- 
tral 7", anti-correlated with the mean velocity which has pinched 
isocontours. The long-slit data of Simien &Prugniellfl997al) con- 
tains a hint of this central kinematical structure. 

NGC 2699: The kinematical major-axis in the central region is 
misaligned from the photometric major-axis. The mean velocity 
map shows a rapidly rotating component with a local maximum 
at about 5" from the centre: within this region, /13 departs signif- 
icantly from zero and is anticorrelated with V. Th ere is a weak 
hint of such a decoupled structure in the data of Simien & Prugniel 
(2000). 

NGC 2768: This galaxy shows a rather cylindrical veloc- 
ity field, with a shallow gradient across the rotation axis 
( Fried & Illingworth 1994). There is a pronounced dip in the 
central dispersion field, also clearly present in the data of 
ISimien & Prugniell jl997b |). Signi ficant dust extinction is present 
north of the center lMicharJl999l) . 

NGC 2974: A very rapidly rotating galaxy which has an /13 field 
with one of the strongest amplit udes. The /14 map exhibits a central 
dip, as found by e.g.. lCinzano & va n der Marel| |l994l). but stays 
positive, in contrast to the data of Be nder et alj <1994). 

NGC 3032: This dusty galaxy shows mild rotation and low dis- 
persion values. 

NGC 3156: Another dusty galaxy with rather low dispersion 
everywhere (a < 75 km s~ ), consistent wit h other authors 
iRix. Carollo & Freemarll999tlBender et alll994l) . 

NGC 3377: The zero velocity curve is slightly misaligned 
with the p hotometr ic minor-ax is of the galaxy (Paper I; 
ICopin. Cretton & Emsellerr]|2004l) . giving the modest minor-axi s 
rotation observed in the long-slit study of lHallidav et alj 1200 ll) . 
and the dispersion map shows a clear elongation along the photo- 
metric minor-axis of the galaxy. 

NGC 3379: This roundish galaxy has /13 values being 
anti-correlated with the mean velocity, and has fairly con- 
stant positive hi values: this is consi s tent with other pub- 
lished data (e.g.! lHallidav et al J 1200 it iGebhardt et ail l200rA 
IStatler & Smecker-Hanell999l) . 

NGC 3384: The dispersion map exhibits a dumbbell-like shape, 
coincident w i th a rapidly r otating component in the central 5" 
lFishedll997l:lBusarello et alJll996l) . with a strong anticorrelation 
signature in /13. The velocities are nearly cylindrical in the outer 
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part of this galaxy classified as SBO (Paper II). Our hi map shows 
a sligh t dip at the centre, consistent with the long-slit profile of 
Fisher il997l) . Our hi map is consistent with hi being positive ev - 
ery where, in contrast to the measurements of Pin knev et alj f2003). 

NGC 3414: A central counter-rotating component is visible in 
the velocity map of this peculiar galaxy (Arp 162), with the h 3 
field showing a corresponding sign reversal just outside this kine- 
matically decoupled structure. The decoupled nature of this central 
component has not been previously reported in the literature, al- 
though long-slit data exist jBertolaet all 19951) . 

NGC 3489: The photometry is perturbed by dust extinction. The 
effect is not evident in the kinematic maps although it may explain 
the local increase in the dispersion map at about 15" West from the 
centre iCaon. Macchetto & Pastoriza 2000). Inside 4" the opening 
angle of the iso-velocities decreases, indicating a fast rotating struc- 
ture. h 3 decreases inwards, with a turn over at a radius of about 4", 
inside of which it shows a small increase. 

NGC 3608: The velocity field displays the prese nce of 

the counter-rotatin g core Ijedrzeiewski & SchechteJ Il988t 
Hallid av et all200ll) in the central 8". As in the case of NGC 3384, 
we measure positive hi values everywher e, in disagreement with 
the data published by Pink nev et al] 12003). 

NGC 4150: The dispersion map exhibits a double peaked struc- 
ture along the major-axis, which corresponds to a region where the 
velocity gradient flattens out. A closer look at the central 5" reveals 
the presence of a counter-rotating structure. 

NGC 4262: This strongly barred galaxy has a remarkably regular 
velocity field, the long-axis of the bar being clearly misaligned with 
the kinematic minor-axis. ^,3 changes sign twice along the major- 
axis. 

NGC 4270: At large radii, the velocity field is close to cylindri- 
cal. Inside 10" the opening angle of the iso-velocity contours de- 
creases, with a corresponding increase in the h 3 amplitude. Th e 
velocity curve in the central ~5 is flat iSimien & Prugnielll997cJ) . 

NGC 4278: This galaxy shows fairly regular kinematics, with 
the mean veloc ity decreasing at radii larger than 10" alo ng the 
major-axis iSchechter & Gunr]|l979tlDavies & Birkinshawl 19881 
The central dispersion peak is elongated along the major-axis. The 
hi map reveals relatively high positive values everywhere except 
inside 6" where it d rops si gnificantl y, consistent w ith the long- 
slit data published b y lvan der Marel & FranxHl993l) . but not with 
Bend er etalJ<1994h . 

NGC 4374: The velocity map is consistent with V ~ ev- 
erywhere in the SAURON field lDaviesl ll98ll) . except for a very 
weak large-scale gradient and a rotation pattern in the central 3". 
With NGC 4278 and NGC 4486, this is among the galaxies in the 
SAURON E/S0 sample with the highest relative contribution from 
emission lines. 

NGC 4382: The velocity field is nearly cylindrical in the outer 
part. It exhibits a strong twist inwards, with the signature of a de- 
coupled compon ent in the central 2", where there is also a sign 
reversal in h 3 iBender et alll994l) . A clear counter-rotation is vis- 
ible from higher resolution OASIS maps iMcDermid et al. 2004). 
The dispersion shows a remarkably strong depression in the cen- 
tral 8-10", clearly seen in lFisheiHl997l) . There is a ring of positive 
hj val ues at about 8", hinted in the long-slit data of Bend eTet alJ 
(1994). 

NGC 4387: This rather boxy galaxy iPeletier et al]|l990h has a 
central dispersion drop, already noticed bv lHallidav et alJfeOOll) . 

NGC 4458: The KD C is clearly visible in the SAURON velocity 
map of this galaxy ("see lHallidav et alT2 001). and the velocity out- 
side 5" is consistent with being equal to zero. 



NGC 4459: The small opening angle of the iso-velocities in the 
central 10" corresponds to a region whe re h 3 is anti-cor related with 
V, and has a rather high amplitude (see IPetersonl 1 978). 

NGC 4473: The velocity dispersion field exhibits a complex 
morphology, with a region of high dispersion along the major- axis 
which w idens at larger radii, and a central drop inside 3", not ob- 
served in B ender et alJll994l) . The widening in the dispersion map 
is particularly significant outside 10" where the mean velocity is 
observed to decrease outwards. This also corresponds to a region 
of lower hi. 

NGC 4477: The SAURON field of view zooms in on the bar of 
this galaxy where the kinematic and photometric major-axis are 
misaligned and the velocity field shows a significant twist of its zero 
velocity curve within the ce ntral 4" and nearly cylindrical rotation 
outside. Ijarvis etal] ll988l) noticed the existence of considerable 
minor-axis rotation as well as a flat velocity dispersion profile with 
a central dip which is not visible in our data. 

NGC 4486: This well-known giant elliptical galaxy has veloci- 
ties consistent with zero everywhere in the SAURON field of view, 
confirming the absence of rotation, and high central velocity disper- 
sion fcarge nt et alll9 7HlDavies & Birkinshaw 1988; Bend er et all 
Il994llvan der Marel 1994 . 

NGC 4526: This SO object has a prominent dust disc which is 
not only visible in the reconstructed SAURON image, but also in- 
fluences the observed kinematics. The dumbbell shape in the dis- 
persion map is due to a fast rotating stellar component, well visible 
in the velocity map. The data of Pelle grini. Held & Ciottil 1 19971) 
also shows the drop in velocity between 10" and 20", which corre- 
sponds to this component. h 3 reverses sign at about 20" along the 
major-axis. 

NGC 4546: The velocity field of this barred galaxy displays a 
regular rotation pattern iBettoni & Gallettal 1 1 997h . anticorrelated 
with h 3 . 

NGC 4550: This object has two counter-rotating stellar disks of 
similar mass 1 Ri x et alJll992f) which, at the spectral resolution of 
SAURON, produce a large region with a near zero mean veloc- 
ity V and a dispersion increasing outwards along the major-axis. 
Note the extraordinary decoupling revealed by SAURON, the outer 
disk being counter-rotating with respect to the main body of the 
galax y further away from the equatorial plane, also noticed by 
lAfanasiev & Sirchenkol 120021) . hi is significantly lower outside 
5" along the major-axis. 

NGC 4552: We detect weak rotat ion in this galaxy a nd a corre- 
sponding weak correlated /13 field lBenderetal"]ll994l) . The zero 
velocity curve is significantly twisted. 

NGC 4564: This flattened galaxy shows the signature of a disk- 
like component in the pinched V map , anti-correlated with h 3 
iBender et al Jl 1 9941 : lHallidav et alj|200ll) . The large h 3 amplitude 
along the major-axi s outside ~ 5" is however not observed by 
|PuiknevetaiT f2003) who also measured significantly negative hi. 

NGC 4570: The inner part of the velocity field and h 3 
map suggest the presence of a fast rotating disk-like struc- 
ture. This c omponen t was also seen in lo ng-slit observations 
fvan den Bosch. Jaffe & van dg^MareJ ^98j ) and interpreted as 
due to bar-driven evolution Ivan den Bosch & Emselleml 19981) . 

NGC 4621: This object is another case showing indication 
for a disk-like compo nent with strong an t icorrel ation between 
V and h 3 IB ended Il990t IBender et al] Il994l) . A counter- 
rotating component is also detected inside 2", previously dis- 
covered wit h adaptiv e optics ass isted integral-field spectroscopy 
IWernli. Emsellem & Copinl2002l) . 

NGC 4660: The velocity map shows the superposition of two 
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rapidly rotating components: an inner one ins ide 5" and a thicker 
structure outside 7". This galaxy was used by B ender et alj 11994^ 
as a prototype for discy ellipticals, due to its pointed isophotal 
shape and the strong hz anticorrelated with V . 

NGC 5198: The velocity field reveals a central KDC, rotat- 
ing nearly perpendicularly with respect to the outer body, both 
being misaligned with respect to the photometric axes of the 
galax y. A "suspect" central v elocity profile was noticed by 
Ijedrzeiewski & SchechteJ Il989h . 

NGC 5308: This is a typical case of a close to edge-on disk 
galax y showing a rapidly rotating component ( Sei fert & Scorzal 
Il99d) and a double sign reversal in hz . The central dispersion struc- 
ture is box-shaped. 

NGC 5813: A galaxy with a well-known KDC 
lEfstafhiou Ellis & CarteJ ll98Ct Il982t iKormendvl Il984t 
lBenderetalJll994 evident in the SAURON velocity map (see 
Paper II). The galaxy has negligible rotation at larger radii, and 
shows weak evidence for minor-axis rotation. The dispersion field 
also shows an unusual ring-like depression at about 5" from the 
center. 

NGC 5831: The SAURON maps needed significant rebinning to 
reach the re quired S/N, but the velocity field rev eals the well- 
known KDC i lDavies et alll983tlPeletier et alJl990I) . 

NGC 5838: This boxy galaxy is a strong rotator, with a dou- 
ble sig n reversal in hz- Minor-axis spectra of Fa lcon-Barroso et alj 
(2003) show a small amount of rotation, very probably due to the 
extinction produced by a dust lane most prominent on the South- 
West side within 3" of the centre. 

NGC 5845: The velocity field of this galaxy shows two distinct 
peaks on each side, at about 2" and 8" from the cente r. These fea- 
tures a re barely detected in the long-slit data of Simie n & Prugniel 
(2002). The central kinematically decoupled structure corresponds 
to a lower dispersion and a peak in hz. 

NGC 5846: This is an example of a giant elliptical with a very 
low (but present) rotation ICarollo. Danziger & Busonlll993l) . and 
a high central dispersion. The SAURON field of view includes a 
foreground star and a companion (NGC 5846A) North and South 
of its nucleus, res pectively, hi i s high everywhere, confirming the 
trend observed bv llBender et alj C?94). 

NGC 5982: The S AURON data confirm the pre s ence of a 
KDC first detected bv lWagner. Bender & Moellenhofj Il988h . The 
SAURON maps reveal that the central rotation axis seems to be 
nearly perpendicular to the rotation of the main body, the central 
hz structure following this orientation. 

NGC 7332: Another strongly boxy galaxy with a KDC (c ounter- 
rotating) in the central 3" (see Falcon-Bar roso et alj |20041. and a 
high central dispersion region elongated along th e major-axis (see 
also data bv [ F isher. Illingworth & Franxl 11994: ISeifert & Scorzj 
ll996tlSimien & Prugniell997bl : ISil'chenkoll999l) . 

NGC 7457: The velocity map exhibits a nearly cylindrical rota- 
tion, as well as a central co unter-rotating structure (see also data by 
ISimien & Prugnielll997bl) . 

APPENDIX B: TESTS OF THE KINEMATICS 
EXTRACTION 

Bl Comparison of pPXF with FCQ 

In Fig.lBTl we show a comparison between the stellar kinematics 
recovered via FCQ and our pPXF technique (see Sect. l2.5.TV We 
illustrate this with the case of NGC 4473 which is free of (detected) 
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Figure Bl. Comparison between the stellar kinematics of NGC 4473 ex- 
tracted via the FCQ (crosses) and the pPXF (solid lines) methods. The major 
and minor-axis profiles are respectively shown in the left and right panels, 
with, from top to bottom V, a, h z and /14. The long-slit major-axis mea- 
surements of Be nder et alj ll994l) are shown (filled dots) for comparison. 

emission lines (see Fig. lB2> . In Fig. lBll we present the major (left 
panels) and minor-axis (right panels) profiles for the first four ex- 
tracted parameters, namely V, a, hz and hi. The agreement be- 
tween the FCQ (crosses with error bars) and pPXF (solid lines) 
methods is excellent except for a slight discrepancy observed in the 
hi profiles. FCQ tends to provide hi values systematically lower 
than pPXF. Although this offset depends significantly on the con- 
tinuum subtraction procedure (e.g., degree of polynomial) and the 
spectral domain (including the region around Hf3 or not) used in 
FCQ, we could not fully reconcile the FCQ and pPXF hi mea- 
surements. This difference in the hi measurements may be related 
to the negative bias in the determination of hi via FCQ, whe n the 
LOSVD is not well sampled, noticed by I Joseph et al] fcOOlh (see 
their Fig. 13) and in lPinknev et alj J2003I) . 

B2 Examples of resulting fits 

Here we present a few illustrative examples of spectral fits ob- 
tained using the pPXF and optimal fitting procedure described in 
Sect. 12. 5l Figure lB2l displays the results obtained for five galaxies 
in the SAURON E/S0 sample, probing a range in absolute luminos- 
ity, absorption-line depths and degree of contamination by emission 
lines: 

• NGC 4150 which has rather deep H/3 absorption lines and 
some contribution from emission lines, 

• NGC 4278 which shows very significant emission line con- 
tamination, 

• NGC 4473, an intermediate luminosity elliptical with no de- 
tected emission lines, 

• NGC 5982 as the brightest galaxy in the SAURON E/S0 sam- 
ple, and 
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Figure B2. Resulting spectral fits obtained with the pPXF routine for a few 
galaxies in the SAURON E/S0 sample. Each panel corresponds to a different 
galaxy (see name label at the left), with the upper spectrum taken at the 
center and the lower one at a radius of about 10". The SAURON spectra 
are drawn as thick grey lines, and the corresponding fit obtained from the 
pPXF routine drawn as thin solid lines. All spectra have been shifted to rest 
wavelengths, as well as normalized and vertically offset for better legibility. 
The expected locations of a few emission lines are indicated at the top of 
the Figure. 



• NGC 7457 for a case of a low dispersion galaxy with detected 
emission. 

In each case, we present a central spectrum as well as a spec- 
trum at a radius of about 10" (chosen to lie at 45 "from the photo- 
metric major-axis). As mentioned in Sect. 12.51 all fits are excellent, 
with local discrepancies between the fitted and observed spectra 
emphasizing the presence of emission lines (H/3, [O III], [N I]). The 
use of an individual optimal template spectrum for each SAURON 
spectrum is crucial as line depths can vary significantly within the 
SAURON field of view. These plots demonstrate the flexibility of 
our procedure. 



^..l-UAyjit - 



iiffwpiLpjm#n^f + + 



0.2 0.4 0.6 

RMS scatter (%) 



-40 L-l. 




0.2 0.4 0.6 0.6 

RMS scatter (%) 



0.30 
0.15 
' 0.10 



0.00 

o.o 

















+ + 





1.2 0.4 0. 

RMS scatter (%) 




0.2 0.4 0.6 o.a 

RMS scatter (%) 



Figure B3. Effect of template mismatch on retrieved stellar kinematics. A 
random perturbation was added to a template spectrum before convolving 
it with a LOSVD parametrized as a fourth order Gauss-Hermite series. The 
input LOSVD parameters were Vi n = 20 kms , a"i n = 200 kras -1 , 
^3 in = 0-1 an d h4 ln = —0.1. The measured kinematics (crosses) for 
1000 different Monte Carlo realizations of the perturbation, are shown as 
a function of the corresponding RMS (in percentage) of the residual mis- 
match between the input spectrum and the best fitting model (see text). The 
thick gray line indicate the true input parameters of the LOSVD. 



sured stellar kinematics as we perturb the spectral shape of either 
the stellar template or the galaxy spectrum. 

We performed this in a way similar to the procedure described 
bv lRix & White! Jl992h . We first chose an optimal template repre- 
sentative of the ones obtained for the galaxies in the E/S0 SAURON 
sample. We added a random Gaussian perturbation to that tem- 
plate before convolving it with a LOSVD parametrized as a Gauss- 
Hermite s eries to simulate a 'g alaxy' spectrum. In contrast to what 
is done in lRix & White 1 1992), we preferred not to add any addi- 
tional Poissonian noise, to isolate the sole effect of template mis- 
match. We then retrieved the kinematic parameters with our pPXF 
routine (see Sect. l2.5.Tl to finally compare the extracted and input 
values. This is repeated for many different Monte Carlo realizations 
of the perturbation. 

In Figure lB3l we show such a comparison for an input ve- 
locity of Vin = 20 km s _1 , a dispersion of a- m = 200 km s _1 , 
/i3,in = 0.1 and h^in = —0.1. The plots show the retrieved kine- 
matics versus the RMS of the residual mismatch between the sim- 
ulated 'galaxy' spectrum and the best fitting model. The maximum 
uncertainty in the retrieved kinematic parameters increases almost 
linearly with the scatter in the fit. In the case of the SAURON data 
presented in this paper, we measured (on the highest S/N spec- 
tra where template mismatch dominates over Poissonian noise), a 
typical residual mismatch in the fit of about 0.4% RMS , with a 
maximum value of 0.55%. The above test then indicates an upper 
limit to any systematic error on I13 and hi of < 0.03, although in 
practice our actual error is likely to be significantly lower than this 
maximum. 



B3 Template mismatch 

The effect of template mismatch is difficult to quantify, as galaxy 
spectra are a complex (luminosity-weighted) mix of different stellar 
populations with different velocity distributions. We can however 
perform some simple tests by examining the change in the mea- 



